A characterization of the diffuse Galactic emissions 
at large angular scales using the Tenerife data 
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Abstract — The Anomalous microwave emission (AME) has 
been proved to be an important component of the Galactic 
diffuse emission in the range from 20 to 60 GHz. To discriminate 
between different models of AME low frequency microwave data 
from 10 to 20 GHz are needed. We present here a re-analysis of 
published and un-published Tenerife data from 10 to 33 GHz at 
large angular scales (from 5 to 15 degrees). We cross-correlate 
the Tenerife data to templates of the main Galactic diffuse 
emissions: synchrotron, free-free and thermal dust. We find 
evidence of dust correlated emission in the Tenerife data that 
could be explained as spinning dust grain emission. 

Index Terms — ISM: general - ISM: clouds - Methods: data 
analysis - Cosmology: observations - Submillimeter - Catalogs 

I. Introduction 

The anomalous microwave emission (AME) is an important 
contributor of the Galactic diffuse emission i n the range 
from 20 to 60 GHz. It was first identified by dLeitch et all 
119971) as free- free emission from electron s with temperature, 
T e > 10 6 K. iDraine & Lazarianl d!998al) argued that AME 
may result from electric dipole radiation due to small rotating 
grains, th e so-called spinning dust. M odels of the spinning dust 
emission dDraine & LazarianL Il998bl) show that an emissivity 
spectrum pea king at around 2 0-50 GH z is able to reprod u ce the 



1997), 



1998), 



observations dFinkbeineri 120031: Ide Oliveira-Cost a et all 
I Watson et all 120051: llglesias-Grothl 120051: ICasassus et all 
120081 : iDickinson et all 120091 : iTibbs et all 120101) . The initial 
spinning dust model has been refined regarding the shape and 
rotational pr o pertie s of the dust grains dAli-Haimoud et all 120091 : 
iHoang et all 12010L 120111 : ISilsbee et all 120111). An alternative 



explan ations of AME was proposed bv IDraine & Lazarianl 
based on magnetic dipole radiation arising from hot 
ferromagnetic grains. Observations have placed limits of a few 
per cent on the fracti o nal polarization t oward s AME t argets 
dBattistelli et all 120061: ICasassus et all 120061: iKogutL 120071 : 
iMason et all 120091 : lLopez-Caraballo et all 120111) . This excludes 
perfectly-aligned single-domain magnetic grains, however other 
alignments and grain compositions prod uce similarly low levels 
of polarization dDraine & HenslevL 120121) . 



A correlation between microwave and in frared maps, mainly 
dominated by dust thermal emission dDesert et al I IT990T) . 
was observed for variou s exper iments, for exam ple on 
COBE/DMR dKogut et all ll996allbT) . OVRO (Leitch et all 



119971) . Saskatoon dde Oliveira-Costa et~ail 
survey at 19 GHz dde Oliveira-Costa et all 
Tenerife dde Oliveira-Costa et all 119991) . A similar signal 
was find in compact regions by dFinkbeinerl [2003) 
and in some molecular clouds b a sed o n data from 
COS MOSOMAS dGenova-Santos et all HoTll: [Watson et all 
120051) . AMI dAmi-Consortium: Scaife et all l2009aUb[) . 
CBI dCasassus et all 1200(3 : ICastellanos et all 120111). 
VSA dTibbs et all 120101) and Planck dPlanck-Collaborationl 
120111) . A recent stud y of the Small M agellanic Cloud also claims 
a detection of AME dBot et all I2010T) . 



Independently, I Bennett et alj (2003) proposed an alternative 
explanation of AME based on flat-spectrum synchrotron 
emission associated to star-forming regions to explain part of 
the WMAP first-year observations. T his hypothesis seems to 



be in disagreemen t with results fr om | de Oliveira-Costa et al 
2004) : iFernandez-Cerezo et alj i20m : iHildebrandt et al 



' 2007l) : lYsard et alJ d2010fT which showed that spinning dust best 



explai ned the excess below 20 GHz. Furthermore, I Dav ies et alj 
12006) showed the existence of important correlation between 
microwave and infrar ed emission in regio ns outside star-forming 
areas. More recently, I Kogut et alJ 1201 ll) discussed the fact that 
spinning dust fits better to ARCADE data (3.8 and 10 GHz) 
than a flat-spectrum synchrotron. 

To discriminate between the different AME models and 
from alternative explanations such as those discussed above low 
frequency microwave data in the range from 10 to 20 GHz at 
different angular resolutions are needed. Indeed in this frequency 
range we expect the AME spectrum to be significantly distinct 
from magnetic dust and flat-spectrum synchrotron. At this 
respect the Tenerife data set, from 10 to 33 GHz and at large 
angular scales (from 5 to 15 degrees), is unique. We present in 
the following a re-analysis of these d ata including previously 
published data I Gutierrez et alj (2000) and un-published data 
since January 1998 to December 2000. The paper is structured 
as follows. Section JI] presents the Tenerife data and discuss the 
re-processing of these data. In Section 1 1 1 M I we describe the main 
Galactic emission mechanisms and the associated templates used 
in the analysis. Section[lV]discusses the point-source contribution 
to the Tenerife data. In Section [V] we present the cross correla- 
tion between the Tenerife data and the Galactic templates are 
presented. The main results are disc ussed in Section [vTl Finally 
we draw conclusions in Section Ivnl 

II. The Tenerife data 

The Tenerife experiment observed unidimensional scans at 
constant declination at 10, 15 and 33 GHz using a triple beam 
pattern of FWHM 5° and of 8.1° beam spacing. The region of 
the sky centered at declination 40° was chosen for observations 
because it corresponds to the largest area of the sky at high 
latitudes where contamination from foregrounds is a minimum. 
To reconstruct 2D maps of the sky, consecutive declinations 
separated by half the beamwidth (2.5°) are observed. Each single 
declination is repeatedly observed until sufficient sensitivity and 
full RA coverage is achieved. The scheduling of the observations 
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Fig. 1. 10 GHz stacked scans. From top to bottom we plot the stacked 
scans for the 10 GHz radiometer from declination 30° to 45°. The data have 
been binned into 4° pixels. 



Fig. 2. 15 GHz stacked scans. From top to bottom we plot the 15 GHz 
stacked scans in the declination range 25° to 47°. 5. The data have been 
binned into 4° pixels. 
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Fig. 3. 33 GHz stacked scans. From top to bottom we plot the stacked scans 
for the 33 GHz radiometer in the triple beam configuration for declinations 
40° to 42.5°. The data have been smoothed into 4° pixels. 



takes into account the position of the Sun so that its contribution 
to the data is a minimum. This requires observations of the 
same declination at different times of the year for full RA 
coverage. Observations were performed day and night. Day-time 
observations present an increase in total power and noise with 
respect to night-time observations due to receiver gain changes. 
In extreme cases, data observed at day time have to be removed. 

The atmospheric contribution to the data depends on 
frequency, being severe at 33 GHz and relatively small at 
10 GHz. The observing efficiency depends mainly on the 
weather conditions. Data strongly affected by atmospheric 
contamination cannot be used. In addition, technical problems 
such as power cutoff, warmed and oscillating HEMT amplifiers, 
RF interference, failure in the electronic systems, etc. can also 
affect the data. At 10 and 15 GHz more than 80% of the 
observed data are useful. However, at 33 GHz only about 10% 
of data are kept, due mainly to atmospheric effects. 

We use in the following the full Tenerife data set which 
includes new data (taken from January 1998 to December 200) 
with respect to previ o usly published release s (see for example 
IGutierrez etail <2000l) : lHancock et alj fl99l) ). These data have 
been completely reprocessed including hand editing to remove 
clearly systematic contaminated regions and accurate removal of 
the atmospheric emission using a MEM based b aseline removal 
technique as described in IGutierrez et alj fcOOOf). The latter has 
been improved to enlarge the sky region for which atmospheric 
residuals are negligible. In particular, this allows us to consider in 
the following analysis regions at low Galactic latitudes. Figures [lj 
H] and |3] present the cleaned Tenerife data at 10, 15 and 33 GHz 
respectively. The main properties of these clean data: central 
declination of the scan, mean temperature, mean noise per beam 
and r.m.s per 1° pixel in the RA interval 150° to 250° are 
presented in Tables JUl] and |m] for the 10, 15 and 33 GHz 
channels respectively. 

III. Galactic diffuse emission. 

We present in this section the main known Galactic diffuse 
mechanisms: synchrotron, free-free, and vibrational and rota- 
tional dust. 

A. Synchrotron 

Synchrotron emission results from cosmic-ray electrons ac- 
celerated in magnetic fields, and thus, depends on the energy 



TABLE I 

10 GHz data. Mean temperature, mean noise per beam and r.m.s 
per 1 ° pixel for the 1 GHz FINAL stacks in the RA INTERVAL 150° 

TO 250°. 



Declination (deg) Mean T (fi K) cr(p K) r.m.s (/.tK) 



30.0 


-3.7 


52.9 


150.5 


32.5 


-5.0 


45.9 


147.8 


35.0 


10.0 


47.1 


140.5 


37.5 


7.5 


51.0 


153.1 


40.0 


13.6 


49.6 


147.4 


42.5 


7.9 


54.1 


146.0 


45.0 


-7.8 


59.2 


215.5 



TABLE II 

15 GHz data. Mean, mean noise per beam and r.m.s per 1° pixel 

FOR THE 15 GHZ FINAL STACKS IN THE RA INTERVAL 150° TO 250°. 



Declination Mean T (fi K) cr(fi K) r.m.s QiK) 



25.0 


3.6 


18.1 


51.7 


27.5 


4.6 


16.3 


48.7 


30.0 


3.4 


19.3 


61.1 


32.5 


3.9 


22.7 


67.2 


35.0 


0.25 


20.6 


74.6 


37.5 


8.3 


22.4 


78.6 


40.0 


7.S 


20.7 


75.7 


42.5 


5.0 


22.6 


81.3 


45.0 


-8.2 


21.3 


75.8 


47.5 


-5.9 


19.3 


70.5 



spectrum of the electrons and the intensity o f the magnetic field 
jRibicki & Lightmaii 11979! : lLongairi. 11994) . The local energy 
spectrum of the electrons has been measured to be, for energies 
contributing to the observed radio synchrotron emission, a power 
law to good approxim ation with index fro m about -2.7 to -3.3 
over this energy range dDe Zotti et aIL Il999). Such an increase of 
energy spectrum slope is expected, as the energy loss mechanism 
for electrons increases as the square of the energy. 

Radio surveys at frequ encies less than 2 G Hz are dominated 
by synchrotron emission jLawson et all 119871) . The only all-sky 
survey available at these frequencies is the 408 MHz map 
jHaslam et ail 119811) . This survey has a resolution of 0°.85 and 
was produced using the Parkes 64-m telescope in Australia for 
the southern sky and the Bonn 100m and Jodrell Bank MK1A 
telescopes for the northern sky. The scanning strategy with the 
Bonn telescope was to fix the azimuth at the local meridian and 
scan up and down in elevation at a rate of about 6°/min. This 
technique reduced the atmospheric contribution to the map but 
led to a set of vertical stripes (constant RA) separated by 7°. 
The quoted errors in the temperature scale are of the order of 
10% and ±3 K in the absolute brightness temperature levels. 

In addition there is the 1420 MHz survey dReich & Reichl 
119881) which covers the declination range 19° to 90° and has a 
FWHM of 0°.58. Stripes are also present in this map due to the 
scanning strategy which consisted of azimuthal scans at constant 
elevation. The errors in the temperature scale are of the order 
of 5% and ±0.5 K in the absolute brightness temperature. 

The 408 MHz and 1420 MHz maps have been used to 



TABLE III 

33 GHz data. Mean, mean noise per beam and r.m.s per 1° pixel 

FOR THE 33 GHZ FINAL STACKS IN THE RA INTERVAL 150° TO 250°. 



Declination 


Mean T (/iK) 


<T(/iK) 


r.m.s (/j,K) 


40.0 


-3.1 


27.2 


80.8 


42.5 


0.2 


43.7 


167.7 



4 




Fig. 4. From top to bottom. 408 MHz synchrotron map. Resolution: 0°.85 Pixels: 0°.33. WHAM Ha map. This map is courtesy of the WHAM 
collaboration. Resolution: 1° Pixels: 1°. Combined DIRBE+IRAS 100/^m dust map. Resolution: 6 arcmin Pixels: 0°.25. 
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determine the synchrotron spectral index at radio frequencies 
jReich & Reicfl, 119881 : iDavies et allU996l) . Assuming T v oc v' 13 
and after destriping and correction for zero levels, spectral 
indexes, /3, in the range 2.8 to 3.2 were found. The spatial angular 
power spectrum of the synchrotron emission C a f vnc is poorly un- 
derstood but is believed to be C^ ync oc t~ z dBouchet & Gisperll 
119991) . At high Galactic latitudes in the region observed by 
the Tenerife experiment the synchrotro n spatial angula r power 
spectrum is slightly flatter than Z 2 dLasenbvl 119961) . In the 
following s ections, we will use the destriped version of the 408 
MHz map jDavies et all 119961) as a template of the synchrotron 
emission. This map is shown on Figure 0] 



B. Free-Free 

When a charged particle is accelerated in a Coulomb field 
it will emit radiation which is called braking radiation or 
Bremsstrahlung. The Galactic free-free emission is the thermal 
bremsstrahlung from hot electrons ( ~ 10 4 K ) produced in 
the i nterstellar gas by the UV radiation field dDe Zotti et all 
119991) . This emission is not easily identified at radio frequencies, 
except near the Galactic plane. At higher latitudes it must 
be separated from synchrotron emission by virtue of their 
different spectral indices, since the spectral index of optically 
thin free-free emission is &'' = 2.1. 

The diffuse Galactic recombination radiation, H a is a good 
tracer of free-free emission since both are emitted in the same 
Warm Ionised Medium (WIM) and both have intensities propor- 
tional to the emission measure, 



EM = I n e n p de 



(1) 



The ratio of free-free brightness tempe rature to the H 
surface brightness I a in R (Rayleigh) is dValls-GabaudL 119981) 



T? f {mK] 
UR] 



10 Av' 



2.14 j 1 0.527-^q0.029/T 4 



(1 + 0.08) (2) 



where v is the observing frequency in GHz, T4 is the temperature 
of the electrons in units of 10 4 K , and the last factor 0.08 
corresponds to helium which is assumed completely ionised and 
creates free-free emission like hydrogen but does not emit H Q 
light. 

Recently a full sky survey of U a light has been released by the 
WHA M (Wisconsin H Q Mapper) collaboration dHaffner et all 
119991) . The WHAM instrument consists of a 6 inch dual-etalon 
Fabry-Perot spectrometer with a narrow band filter of FWHM 
~ 20 Awhich images onto a cryogenically cooled 1024 x 1024 
CCD. The resolution of the survey is ~ 1°. The spatial power 
spectrum of free-free emission C{ has not yet been determined 
from the WHAM map howe ver analysis of II , images of 
the North Celestial Pole Area jVeeraraghavan & DaviesL [1997) 
suggests G{* oc r 2/27±oxn . 

In Figure [4] we present a map of the northern sky produced by 
the WHAM survey. In the following sections we will use this map 
as a template for the free-free emission. At intermediate Galactic 
latitudes (say \b\ > 30°) about 10 % of the H Q light is absorbed 
by dust and therefore estimates of free-free emission from H a will 
be systematically lower than the true value. At latitudes below 
10° this correction becomes increasingly uncertain. 



C. Vibrational Dust 

At the higher frequency range ( > 100 GHz) of the microwave 
background experiments, dust emission starts to become domi- 
nant. Dust grains are heated by interstellar radiation, absorbing 
optical and UV photons and emitting energy in the far infrared. 



The intensity of the emission from an ensemble of dust grains 
is given by 

= J <v)di (3) 

where e(y) is the emissivity at frequency v, and the integral 
is along the line of sight. In the Rayleigh-Jeans regime and 
assuming a constant line-of-sight density of dust, 



(4) 



where T b is the brightness temperature of the dust emission. 

The spectrum of the dust emission has been measured at 
millimeter and submillimeter wavelengths by the Far-Infrared 
Absolute Spectrophotometer (FIRAS) and can be fitted by a 
single greybody spectrum of t emperature 17.5 K and emissivity 
oc v 2 jBoulanger et all [j.996) at hi gh Galactic la t itudes . From 
IRAS observations of dust emission iGautier et all 119921) . it was 
found that the spatial power spectrum of the dust fluctuations 
is Cf U3t oc £~ 3 . This has also been confirmed at larger angular 
scales by the COBE-DIRBE satellite. 

In the following secti ons we will use the combined IRAS- 
DIRBE map at 100 ^m dSchlegel et all 119981) as a template for 
the dust emission. This map has resolution of FWHM ~ 6 arcmin 
and covers the full sky. Zodiacal light and point sources have been 
removed from the map and the regions of the sky which were 
not observed by the IRAS satellite have been replaced by DIRBE 
data. The combined map preserves the DIRBE zero point and 
calibration. This map is in units of MJy/sr. In Figure[4]we present 
the northern sky part of the combined IRAS-DIRBE map. 

D. Rotational Dust. 

Small spinning interstellar dust grains containing 10 2 -10 3 
atoms can produce detectable rotational emission in the 10- 
100 GHz range. This emission depends on the component of 
the electric dipole moment perpendicular to the angular velocity 
of the gr ain and on the phys i cal p roperties of the interstellar 
medium (Drajne & Lazarianl 119981) . For these small grains, 
rotational excitation is dominated by direct collisions with ions 
and plasma drag. The very smallest grains (N < 150) have their 
rotation damped primarily by electric dipole emission; for 150 
< N < 10 3 plasma drag dominates. 

In the following sections, we show models for the spectrum of 
the spinning dust emission provided by Dr. Draine (private com- 
munication). These models depend on a large set of parameters 
such as the distribution of grain sizes, the charge of the grains, 
the composition of the grains and the physical properties of the 
interstellar medium which were fixed by the authors. However, 
the normalization of the model can be assumed a free parameter 
although it is related to the hydrogen column density in the 
interstellar medium which is of the order of a few times 10 20 . 

The spatial distribution of small dust grains is not well-known 
although it seems reasonable to believe that it is not different 
from that of larger grains but for dense regions where dust 
coagulation may deplete small grains. For the purpose of this 
work we will use the IRAS-DIRBE 100 /im map as a template 
for the spinning dust grains. 

IV. Galactic and Extra-Galactic point sources. 

The contribution from resolved point sources to the Tenerife 
data at 10 and 15 GHz was extrapolated from data of the 
Michigan monitoring program (M. Aller and H. Aller 2000, 
private communication). The Michigan program regularly mon- 
itors point sources with fluxes above 0.5 Jy at 4.8, 8.0 and 
14.5 GHz. The Michigan catalog is neither complete in flux or 
time domain. Well-known strong variable sources such as 3C345 
are regularly observed and data at all observing frequencies 
are available. Weak sources are poorly observed and often 
data are available at a single frequency of the three possible. 
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TABLE IV 

Contribution from unresolved point sources to the Tenerife 
experiment. 



v (GHz) 




10 


< 7.5pK 


15 


< 4.0/iK 


33 


< 0.6/iK 



Further, data are available up to June 1999 while the Tenerife 
experiment operated until September 2000. For the 33 GHz 
Tenerife data we used the Metsahovi catalog (Metsahovi group, 
private communication) which regularly monitors sources above 
1 Jy at 22 and 37 GHz. This catalog is neither complete in flux 
or time although observations are available up to January 2000. 

We have developed software to produce time and frequency 
complete point source catalogs. The flux of the sources for each 
frequency was interpolated in time by fitting Fourier Series to the 
data. If the number of independent observations per frequency 
was smaller than 10, the sources were assumed not variable and 
the data was fitted to a constant with time. The extrapolation of 
the source fluxes into the future time was performed assuming 
no variability since last observed data point. Spectral indices 
were calculated for each source from the available data. If only 
observations at a single frequency were available we assumed a 
flat spectrum. We constructed four point source catalogs at 10, 
13, 15 and 33 GHz covering the time range January 1984 to 
December 2000 with one Julian day sampling and a flux limit of 
about 0.5 Jy. 

For each single day of observations with the Tenerife radiome- 
ters we estimated the total contribution from point sources to 
the data. This contribution was calculated from the extrapolated 
catalog at 10, 15 and 33 GHz a nd the flux, S, was conve rted into 
antenna temperature, T, using jRohlfs & WilsonL 120001) 



where fl is the beam solid angle and A the wavelength. The daily 
point source contributions were stacked in the same manner as 
the Tenerife data and subtracted from the Tenerife stacks. In 
the left column of figure [5] we plot the Tenerife stacked data 
at 15 GHz for the 10 declinations presented in this paper and 
overplot in green the contribution from point sources to the data. 
In the right column, we plot the stacks after subtraction of the 
point sources. 

1) Unresolved sources.: A study of the contribution of 
unresolved point sources (ie. weak point sources not detected 
individu ally) to CMB experiments has been produced by 
iFranceschini et aii d!989h . They used numerous surveys, includ- 
ing VLA and IRAS data, to put limits on the contribution to 
single beam CMB experiments by a random distribution of point 
sources. The contribution from unresolved point sources to the 
Tenerife data deduced from the previous analysis is presented in 
table IIVI We consider sources with flux under 1 Jy at 5 GHz 
and a main beam of 5°. 

V. Assessment of the Galactic contribution to 
the Tenerife data. 

To assess Galactic contribution to the Tenerife data, we 
correlated observations with the Galactic synchrotron, dust and 
H a emission templates described in the previous sections. We 
convolved the foreground maps with the Tenerife beam at each 
of the Tenerife frequencies before the correlation was performed. 
The Tenerife data used in the correlation are the final stacked 
data presented in Section |Hj At 15 GHz we used ten declination 



stacks in the range 25° 47°. 5; seven declination stacks at 
10 GHz in the range 30° — 45° ; and only two declination stacks 
at 33 GHz covering declinations 40° and 42°. 5. The data were 
processed so that data at |6| < 15 were excluded from the baseline 
fit and therefore were neither stacked or reconstructed. The 
discrete point sources were subtracted from the Tenerife stacks 
as discussed in the previous section. 

A. The Method 

To simultaneously correlate the Tenerife data to the three 
Galactic te mplates we use a met hod which was first applied to this 
problem bv lGorski et alj d!996T> to fit Galactic and extra-Galactic 
templates in Fourier space to the COBE-DMR da ta. This me thod 
was a pplied to the Tenerife data first by | de Oliveira-Costa et alj 
H999I) and then bv lMukherjee et alJ fcOOll) to study the possible 
emission of spinning dust at the Tenerife frequencies. 

Assuming that the microwave data is a superposition of CMB, 
noise and Galactic components, we can write 

y = aX + xcmb + n (6) 

where y is a Tenerife data vector of N pixels; X is an N x M el- 
ement matrix containing M = 3 foreground templates convolved 
with the Tenerife beam; a is a vector of size M that represents 
the levels at which these foreground templates are present in 
the Tenerife data - correlation coefficients for each foreground 
template; n is the instrumental noise in the data; and xcmb is 
the CMB component of the data. For this analysis we assume 
the noise and CMB to be uncorrelated. 

The minimum variance estimate of a is given by 

a = [x T C* _1 x] 1 X T C' 1 y (7) 
with errors given by as, t = J^i/ 2 where Yl i s given as 

£ = <a 2 )-(a> 2 = [x^-'xy 1 (8) 

In the above, C is the total covariance matrix, the sum of the 
noise covariance matrix and CMB covariance matrix. The noise 
covariance matrix of the Tenerife data is taken to be diagonal 
- no correlation between different pixels. The CMB covarian ce 
matrix was obtained analytically following Zaldarriagaj U995t) . 

This correlation method produces minimum variance and 
unbiased estimates of a if the following holds 

• The error in the Tenerife data is Gaussian and with zero mean. 
To a very good approximation the instrumental noise in the 
Tenerife data is uncorrelated. 

• The templates perfectly trace the foreground emissions at mi- 
crowave frequencies. There may be components of emission 
present in the Tenerife data, apart from those we have 
identified, for which there are not obvious counterparts 
at other frequencies. Moreover, the templates we use do 
not perfectly trace the microwave emission of the Galactic 
foreground. The latter is specially important in the case 
of vibrational and rotational dust for which the template 
comes from data at much higher frequencies. For example, 
iDraine & Lazarian j!998h proposed a 12/im map would be 
a much better tracer of rotational dust emission because 
small grains will also emit at this frequency. Also, the H Q 
emission is absorbed at low Galactic latitudes by interstellar 
dust and therefore does not fully trace the free-free emission 
at those latitudes. Moreover, the 408 MHz map has residual 
striping. 

• The templates are perfect - error bars equal to zero. The error 
estimated for the templates are about 5 - 10 % of the signal. 
They are not corrected for in this analysis because they are 
small compared to the errors in the Tenerife data. 

• The templates are not correlated. It there is correlation 
between the different templates the set of minimum variance 
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Fig. 5. Extrapolated point source contribution to the 15 GHz Tenerife data. Left column: observed scans. Right column: point source-corrected scans. 
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TABLE V 

Cross-correlation results. See text for details 



Galactic cut 


V (GHz) 




a ± a- 


AT (/iK) 


b > 20° 


10 GHz 


4QB [\,] | |/ 




17 ± 6 


17.4 ± 6.0 


100 fin 


0.3556 


23 ± 17 


7-5 ± 5.5 


WHAM 


0.551:? 


42 ± 20 


10.5 ± B.5 


15 GHz 


408 MHz 


1.0244 


9 ± 3 


9.2 ± 3.0 


100 fim 


0.3308 


66 ± 8 


21.8 ± 2.6 


WHAM 


0.2758 


— 10 ± 10 


— 2.7 ± 2.7 


33 GHz 


40s MHz 


0.8904 


4 ± y 


3.6 ± 8.0 


100 M m 


0.301$ 


86 ± 27 


26.0 ± 8.0 


WHAM 


0.?76(4 


21 ± 51 


5.8 ± 1U.U 


b > 30° 


10 GHz 


408 Ml 1/ 


— q 8722 — 


15 ± 7 


13.0 ± 6.0 


100 ura 


0.1018 


V2 ± so 


14.6 ± 6.0 


WHAM 


0.1435 


83 ± 42 


11.9 ± 6.0 


15 GHz 


408 Mil/. 


0.8544 


13 ± 4 


11.1 ± 3.4 


100 Lira 


0.2036 


-1 ± 14 


-0.2 ± 2.8 


WI 1AM 


0.1504 


30 ± 20 


4.5 ± 3.0 


33 GHz 


408 MHz 


0.7749 


7 ± 10 


5.4 ± 7.7 


100 nm 


0.2143 


yo ± 45 


19.2 ± 9.6 


WHAM 


0.1473 


17 ± 63 


2.5 ± 9.3 


b > 40° 


10 GHz 


408 MHz 


0.7676 


11 ± 9.5 


8.4 ± 8.0 


100 Lira 


0. 1 898 


110 ± 45 


20.9 ± 8.5 


WHAM 


0.1292 


77 ± 52 


9.9 ± 6.7 


15 GHz 


408 MHz 


0.7415 


15 ±5 


11.1 ± 3.7 


loo 


0.1782 


-11 ± 21 


-1.9 ± 3.7 


WHAM 


0.1354 


51 ± 25 


e.g ± a.4 


33 GHz 


4()8 MHz 


0.6630 


4 ± 15 


a.e ± 9.9 


100 Lira 


0.2138 


66 ± 55 


14.8 ± 11.8 


WI 1AM 


0.1499 


47 ± 6y 


7.0 ± 10.3 



correlation coefficients a is degenerate and therefore we 
would not be able to discriminate among the different fore- 
ground components. At low Galactic latitudes, the Ga lactic 
plane emission dominates and lMukheriee et alj {2001) have 
found the foreground templates are correlated in this region. 
• The correlation coefficient a is the same throughout the area of 
sky for which the correlation is performed. If it is not, the error 
associated with a is systematically underestimated. As a 
double check, we divided the sky observed, into independent 
areas and performed the same correlation test in each, 
calculating a mean correlation coefficient and the dispersion 
of the individual values which in most cases was in good 
agreement with the error of a calculated for the total area. 

B. Cross-Correlation results. 

The cross-correlation results are presented in table |V] The 
correlation was performed for three different Galactic cuts 
b > 20, b > 30 and b > 40. o ga i represents the r.m.s. of 
the Galactic maps after convolution with the Tenerife beam, a 
and £T a are the correlation coefficient and the error associated 
with it. They have units of fiK/K, ^K/(MJy sr ^ and ^K7R 
for the correlation with the 408 MHz, 100 (im and WHAM 
maps respectively. AT is the r.m.s. contribution from the Galactic 
foregrounds to the Tenerife data, which was obtained as o ga i x a. 
This analysis is improved wit h re spect to previous analys es by 
Ide Oliveira-Costa et alj J1999T> and lMukheriee et alj fcoo'll) first, 
because we present Tenerife data at 10 and 15 GHz for a much 
larger area of the sky and with improved sensitivity, and second, 
because we also include in the cross-correlation a template for 
the free-free emission. Moreover, data at 33 GHz have also been 
included in the analysis although the area of the sky observed is 
significantly smaller than at 10 and 15 GHz and consequently the 
error bars in the estimated cross-correlation coefficients much 
larger. We have pictorially summarised the cross-correlation 
results in figure |6j In the top row of the figure we plot the 
r.m.s. contribution from synchrotron to the Tenerife data at the 
three Galactic cuts analysed. In black, we overplot the expected 
synchrotron contribution for a spectral index of —3.0 and derived 
from the r.m.s. level in the 408 MHz map. The contributions are 
systematically lower than expected at 10 GHz and the spectral 
index is much flatter such that for example the contribution at 
15 GHz is consistent with that at 10 GHz for b > 40°. We have 
derived from the Tenerife data spectral indexes of 1.7712, q> 
-0.40±l;g and -0.4O± 2 {% at b > 20°, b > 30° and b > 40° 



respectively, and the r.m.s contributions based on these spectral 
indices is overplotted in red. 

In the middle row we plot the r.m.s. contribution from free- 
free emission to the Tenerife data. In black, we overplot the 
expected free-free emission at microwave frequencies derived 
from equation [2] which are actually in very good agreement with 
the observations. In red, we plot the free-free contribution for 
spectral indexes of -0.5, -2.20 and -1.0 at b > 20°, b > 30° 
and b > 40° respectively. 

In the bottom row of figure [6] we plot the r.m.s. contribution 
from dust emission to the Tenerife data which is significantly 
larger (few orders of magnitude) than expected from vibrational 
dust. In solid black we overplot a r.m.s brightness temperature 
expected from the CNM model of rotational dust proposed by 
iDraine & Lazarianl U998I) . The models were rescaled to fit the 
data. The intensity spectrum for this model peaks at 50 GHz 
which corresponds to a peak in brightness temperature around 
20 GHz of about 300 fiK. From the best fit model to the data we 
have estimated fluctuations in the temperature of the rotational 
dust of ~ 8% at the angular scales of the Tenerife experiment. 

The r.m.s. contributions from dust to the 15 GHz data at 
b > 30° and b > 40°, as well as the contribution from free- 
free at b > 20° seem to be underestimated. This could be 
caused by correlations between the synchrotron, dust and free- 
free templates which would reduce the validity of the minimum 
variance solution and could bias the estimates of the correlation. 
In flgurelll we plot the following correlations, from top to bottom 
synchrotron vs dust; free-free vs dust and free-free vs synchrotron 
for b > 20°, b > 30° and b > 40°. We observe no correlation 
between free-free and synchrotron, moderate correlation between 
free-free and dust and a quite strong correlation between syn- 
chrotron and dust at low Galactic latitudes weakening down at 
high latitudes. Note that at b > 20°, the correlation plots show 
negative vs. negative points which do not follow the correlation 
pattern of the main body of points. The former points correspond 
to the lowest Galactic latitude data which run into a negative 
beam of the triple beam pattern for the Galactic plane crossing. 
The observed correlations between templates can not justify the 
observed lack of dust correlation at 15 GHz data at b > 30° 
and b > 40°. In the following we concentrate on the b > 20° 
region for which we detect significant dust correlated signal at 
all Tenerife frequencies. 

VI. Dust-correlated emission. 

iKogut et all H99d) cross-correlated the COBE Differential 
Microwave Radiometer (DMR) maps with DIRBE far-infrared 
maps and discovered that statistically significant correlation 
did exist at each DMR frequency, which was inconsistent with 
vibrational dust alone. This extra correlation was explained at 
the time as free-free emission. Following this publication, other 
authors have cross-correlated CMB data sets with dust templates 
and found an excess of correlation which has been interpreted 
as free-free emission, flat-spectrum synchrotron or emission from 
spinning dust. In table |Vl] we present an up-to-date list of CMB 
data for which the correlation with dust has been performed. 
In all cases, a multi-template method has been used to perform 
the correlation but no free-free template has been used. We have 
als o added in this table pr e vious anal yses of the Teneri f e data 
by Ide Oliveira-Costa et ail <1999h and iMukheriee et alj fcOOll) 
which did not include the 33 GHz data and covered a much 
smaller area of the sky at 10 and 15 GHz. 

The correlation coefficients we deduced are significantly 

smaller at \b\ > 20 than those calcula ted by 

Ide Oliveira-Costa et alj jl999h and IMukheriee et alj <200l) . 
This is probably due to 50 % larger area (localized signals 
are diluted) and a more careful subtraction of baselines at low 
Galactic latitudes in our data. The inclusion of an extra free-free 
template seems to play minor role on this at 10 GHz and none 



9 




Fig. 6. Spectra of Foregrounds at the Tenerife frequencies. We represent the cross correlation coefficients as a function of frequency for b > 20, b > 30 
and b > 40 (left, middle and right column respectively). From top to bottom, we we plot the r.m.s. contribution from synchrotron, free-free and dust to the 
Tenerife data. In solid black and red we show standard and best-fit models of the electromagnetic spectrum for each Galactic emission. The models for dust 
are rescaled to fit the data. Details on the best-fit models can be found on the text. 
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Fig. 7. Correlation between the foreground templates. From top to bottom, we plot synchrotron versus dust, free-free versus dust and free-free versus 
synchrotron. The correlation is presented for b > 20, b > 30 and b > 40 (left, middle and right column respectively). 
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TABLE VI 

Cross-Correlation results of CMB data with dust. Available cross-correlation coefficients of CMB data sets with dust templates at 

100 fj,M. 



Experiment 


Frequency (GHz) 


\b\ > (degrees) 


a ± a & ^KAMJy.sr L ) 


References 


COBE DMR 


31.5 


20 


18.0 ±2.5 


Kogut et al. (1996) 


COBE DMR 


31.5 


30 


14.5 ± 6.0 


Kogut et al. (1996) 


COBE DMR 


53.0 


20 


6.8 ± 1.4 


Kogut et al. (1996) 


COBE DMR 


53.0 


30 


6.43.4± 


Kogut et al. (1996) 


COBE DMR 


90.0 


20 


2.7 ± 1.6 


Kogut et al. (1996) 


COBE DMR 


90.0 


30 


4.6 ±3.9 


Kogut et al. (1996) 


Saskatoon (Ka band) 


30.0 


NCP 


15.0 ±8.1 


de Oliveira-Costa et al. (1997) 


Saskatoon (Q band) 


40.0 


NCP 


11.8 ±10 


de Oliveira-Costa et al. (1997) 


19 GHz survey 


19.0 


20 


38.5 ± 3.5 


de Oliveira-Costa et al. (1998) 


19 GHz survey 


19.0 


30 


47.1 ± 9.0 


de Oliveira-Costa et al. (1998) 


OVRO 


14.5 


NCP 


209 


Leitch et al. (1997) 


OVRO 


32.0 


NCP 


36 


Leitch et al. (1997) 


PYTHON V 


40.3 




-3.0 ± 18.0 


Coble et al. (1999) 


Tenerife 


10.0 


20 


49.8 ±11.0 


de Oliveira-Costa et al. (1999) 


Tenerife 


10.0 


30 


-8.3 ±31.0 


de Oliveira-Costa et al. (1999) 


Tenerife 


10.0 


40 


84.0 ±54.0 


de Oliveira-Costa et al. (1999) 


Tenerife 


15.0 


20 


71.8 ±4.5 


de Oliveira-Costa et al. (1999) 


Tenerife 


15.0 


30 


94.9 ±15.0 


de Oliveira-Costa et al. (1999) 


Tenerife 


15.0 


40 


72.0 ±33.0 


de Oliveira-Costa et al. (1999) 


Tenerife 


10.0 


20 


71.0 ±18.0 


Mukheriee et al. (2001) 


Tenerife 


10.0 


30 


-7.0 ±32.0 


Mukheriee et al. (2001) 


Tenerife 


10.0 


40 


28.0 ±39.0 


Mukheriee et al. (2001) 


Tenerife 


15.0 


20 


91.0 ±11.0 


Mukheriee et al. (2001) 


Tenerife 


15.0 


30 


29.0 ±20.0 


Mukheriee et al. (2001) 


Tenerife 


15.0 


40 


3.0 ±26.0 


Mukheriee et al. (2001) 


South Pole 94 


30.0 


40 


20.0 ±36.0 


Hamilton & Ganga (2001) 


South Pole 94 


40.0 


40 


68.1 ±42.4 


Hamilton & Ganga (2001) 



at 15 GHz. 

The cross-correlation results for the Tenerife experiment 
presented in this paper confirm the existence of extra dust- 
correlated emission at microwave frequencies and suggest that it 
is not due to free-free emission. The moderate correlation found 
between the dust and free-free templates could confuse the results 
obtained but clearly cannot account for all the observed dust- 
correlated component. A more detailed study is needed to take 
into account correlations between templates in the calculation of 
the correlation coefficients. 

In figure [8] we plot the dust cross-correlation coefficients ( 
|6| > 20) for the Tenerife (diamonds) and COBE-DMR (triangles) 
data. The rest of the experiments presented in table |Vl] were not 
included in this plot because they observe at different angular 
scales and different areas of the sky (see table for details) and 
therefore no direct comparison with the Tenerife and COBE data 
is possible. The Tenerife data suggest a peak in the spectrum 
at about 30 GHz although this is mainly based on the data at 
33 GHz which is significantly much noisier and cover a much 
smaller part of the sky. However DMR data combined with 
the 10 and 15 GHz Tenerife data points prefer a peak in the 
range 15-20 GHz. I n color, we overplot the s ix spinning dust 
models proposed by iDraine & Lazarianl H9981) combined with 
the standard vibrational dust model described in the previous 
sections, for which the amplitude has been taken as a free 
parameter and fitted to the COBE and Tenerife 10 and 15 GHz 
data. We observe that the shape of the spectrum given by the 
data is similar to some of the predicted spectra although the 
models do not indicate a sharp rise in the range 10-15 GHz. 

VII. Conclusions 

We have presented in this paper a re-analysis of the Tenerife 
data including previously published data and new data from Jan- 
uary 1998 to December 2000. This analysis leads to evidences for 
an excess of dust-correlated emission at microwave frequencies 
in the range 10-33 GHz and at large angular scales, from 5 




Fig. 8. Cross-correlation coefficients for the dust-correlated emission. We 

plot the Tenerife data as diamonds and the COBE DMR data as triangles. The 
color solid curves overplotted correspond to joint vibrational and rotational 
dust models as described in the text. 



to 15 degrees. This correlation can not all be associated with 
free-free emission. A combination of the Tenerife and COBE 
data suggests spinning dust emission could account for the extra 
correlation. However, the scatter observed in the data and the 
discrepancy in the spectrum shape indicate other components 
may also be responsible for the extra correlation. Furthermore, 
the analysis does not take into account correlations between the 
Galactic emission templates and this can bias the estimate of 
the correlation coefficients. To correct for this a more detailed 
analysis, which will account for the cross-correlation between 
templates, is needed. 
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